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1. Introduction 

A disorder parameter for dual superconductivity of gauge theory vacuum has been developed 
[|l|] [^][||][^. It is the vev , (jJ.) , of an operator carrying non zero magnetic charge ^. The 
Euchdean version is 

=exp[-^ J d^yE{y,t)b^{x-y)] (1.1) 

Here b± (r) = ^.^"^^ ^ is the field of a monopole in the transverse gauge with 
Vb^=0 

V Ab^ = ^-4nd{n.r)n5^{r^). 

The field E± is the conjugate momentum to the transverse component of the potential A± so 
that jj. is nothing but a translation operator of A± , or 

H{x,t)\A^{y,t)) = \A^{y,t) + ±b^{x-y)) (1.2) 

It just creates a monopole. 

One of the factors ^ at the exponent of Eq(l.l) comes from the Dirac quantization condition 
for magnetic charge, the other one from the fact that the electric field in the lattice formulation has 
an additional multiplicative factor g . The operator /i can be written in the form 

At = exp [-j3A5] (1.3) 

with the usual notation j8 = ^ and AS = J d^yE{y,t)bj^{x — y). As a consequence (jJ.) is the 
ratio of two partition functions and (/i) = 1 at j3 = 0. 

, , Z(S + AS) 

For compact U{1) gauge theory a few theorems have been proved: 

(1) ;U(^,f) is a gauge invariant, Dirac like , magnetically charged operator , and obeys cluster 
property 

(2) (/i) 7^ for /3 < j3f where there is confinement , (/x) = for j8 > j8f i.e. in the deconfined 
phase. j8c is the critical point. 

In ?7(1) gauge theory confinement is therefore produced by condensation of monopoles i.e. by 
dual superconductivity of the vacuum. 

Instead of (jj.) it proves convenient to use the quantity 

(..« 

From Eq(1.4) it follows that p = {S)s — {S + AS)s+as , the subscript indicating the action used in 
the statistical weight. Moreover because of the boundary condition = 1 at j8 = 0, 

(Ai)=exp[/%(/3')^/3'] (1.6) 
Jo 

For SU (N) gauge theories, with and without quarks, N — I operators can be defined (a = 
l,2,...N — 1), and the corresponding order parameters (/x")[0],[|3|]. The definition of jx" has the 
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same form as pL in Eq (1.1) with the field strength replaced by E"^{y,t) = Tr[(^"E{y,t)\\^ where 
<I>'^ selects the direction of the residual abelian gauge field in the abelian projected gauge. For a 
detailed discussion see Ref's[^,[^. The choice of the abelian projection is irrelevant [§][^[|8|]. 

Fig(l) shows the numerical determination of (/x) for compact U{\) gauge theory, Fig(2) shows 
the corresponding quantity p , which presents a strong negative peak at the critical point j3c . The 
analysis goes as follows[^,[^: 
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Figure 1: (/i) for U(l) lattice gauge theory Ref[^. 

1) For p < p" tends to a finite hmit in the thermodynamical limit V = — > oo , and by use 
of Eq(1.6) (At) 

2) For pypcP"^ -\c\L, + c' with c 7^ ,or , by use of Eq(1.6) limi^^ooin) =0 

3) At p « pc the correlation length goes large compared to the lattice spacing and the scaling 
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Figure 2: p versus j3 Ref.[|l| 
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Figure 3: Strong coupling behavior of p at various lattice sizes and am = 0.1335 Ref.|[l 
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Figure 4: Volume dependence of p in the deconfined phase for different values of the magnetic charge.Ref 
[|1 



law holds 



(1.7) 



where T = 1 — j-. v is the critical index of the correlation length , and is typical of the universality 
class of the transition. For weak first order = 3 [ Quenched SC/(3)[|] ,Nf = 2 QCD for 
3d-Ising = 1.6 [ SU(2) [|]], for 3d-0{4) ^ = 1.336. 

The properties 1), 2), 3) as observed in different systems are shown in Figs (3), (4), (5) 



The question we address [13| in this paper is whether (jj.) can be computed in the frame of the 



Stochastic Vacuum Model of QCD[]14|][]15[] [ ]16[ ] . The model consists in expressing physical quan- 
tities in terms of gauge invariant correlators of field strengths , making a cluster expansion of them 
and keeping only the two point cluster. The model provides a good description of many aspects of 
QCD and it would be interesting to know if the distinction between confined and deconfined phase 
could be read in the behavior of the correlators. 
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Figure 5: Scaling of p assuming first order for the deconfining transition. Nf = 2 QCD. Ref. 



2. Cluster expansion of (/i") 



The series expansion of the exponential in Eq(l.l) reads [13] 



(2.1) 



dhi J d\2...J d\nb'l{x-y,) b'l{x-yn){E^^{yut) El{y„,t)) 

According to the stochastic vacuum model one performs a cluster expansion of the correlators: the 
one point cluster is zero by symmetry, and clusters of order higher than 2 are neglected. Keeping 



the correct combinatorics into account [13| the net result is 



(At) = exp[-^ / d'y,d'y2b'l{n -x)b'l{y2 -x){Ek (3^2,0)] 



(2.2) 



Higher clusters are 0{q^) at the exponent. Here E" = Tr [<!>"£]. The gauge invariant correlator at 
the exponent of Eq(2.2) 



{Ei{n,t)Ei{y2,t))=%Jn-y2) 



(2.3) 



in principle depends on the path C used to parallel transport from y\ to y2 but this dependence is 
irrelevant to the study of the ultraviolet and infrared behavior. Since j3 = ^ , by use of Eq's (1.5), 
(2.2) we get 



P" = ^[-^^ Id'y^ I db2bl{n)bl{y2)%i,{n-y2)] 



(2.4) 



In Fourier transform 



(k) = ,H{k) = -4— 

k{k — k.n + ie) k (jik — ie) 

^Uk) = ik^5tj-Kikj)fik^) 



(2.5) 
(2.6) 
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independent on a [16]. 

{k^5ij - kkj)b'^(k)b{{-k) = \H(k)\^ = 4 ~ ^ ^^'^^ 
The equation follows for p 

In the deconfined phase the perturbative expression can be used = ^ and 

pa ^ ^ [-V2L, + 21n(Lv) + constant] (2.9) 

with Ls the spatial size of the lattice , and 

(Lm) IR — cutoff 

L, = = ^ (2.10) 

a UV — cutoff 

In the thermodynamical limit L, ^ oo ,p ^ — oo as in Fig(4), and by use of Eq(1.5), /i ^ 0. 

In the confined phase one expects the same UV behavior, which is dictated by OPE ,but aL^ 
will be replaced by some IR cutoff A so that ,at fixed lattice spacing p" is volume independent as 
in Fig(3) . 

This will never be the case if Eq(2.8) holds, no matter how IR well behaved is the correlator: 
the term oc ^ due to the Dirac string will always diverge. 

This means that the stochastic approach is inadequate in the confined phase. Indeed in that 
phase the vacuum is a Bogolubov-Valatin superposition of states with different magnetic charge 
and the operator n will connect sectors differing by q units of magnetic charge : The Dirac string 
will then end on an antimonopole and the integral will be IR cut-off by a massive propagator and 
be volume independent. 

This can be checked in U{\) theory in the dual formulation of Polyakov [|l7|].The potential 
of the dual field x is there proportional to cos{x) which in the weak coupling is equivalent to a 
mass term, and gives a gaussian distribution like the stochastic vacuum model in QCD. In the 
strong coupling regime the tunneling between the minima of cos{x) provides a Bogolubov-Valatin 
vacuum. 

Like the Polyakov line the order parameter (/i) is singular in the continuum limit a — > 0, but 
its IR behavior at any finite UV cutoff detects confinement or deconfinement. 
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p parameter for different spatial and temporal sizes 

SU(3) pure gauge theory 
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Figure 6: Check of Eq(2.9) 
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